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A lattice of optical ring resonators can exhibit a topological insulator phase, with the direction 
of rotation in the resonators playing the role of spin. This occurs when the inter-resonator coupling 
is sufficiently large, and the synthetic magnetic vector potential set up by the couplers is zero. 
Using the transfer matrix method, we derive the band structure, phase diagram, and the projected 
band diagram showing the existence of spin-polarized edge states. When PT (parity /time-reversal) 
symmetric gain and loss are introduced, the system functions as an optical diode which does not 
require optical nonlinearities. 
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The idea that photonic modes could have non-trivial 
topological properties, analogous to topological states of 
quantum matter, originated some years ago with Hal- 
dane and Raghu [U [2]. These authors argued that, in 
a two-dimensional (2D) photonic crystal with broken 
time-reversal symmetry, there should be modes analo- 
gous to the electronic modes of a "zero-field" quantum 
Hall gas [3|. This was soon confirmed experimentally, 
using a microwave-scale gyromagnetic photonic crystal 
[U [5] . The most striking feature of this system is the ex- 
istence of topologically-protected one-way photonic edge 
states, which can be used for optical isolation; unfortu- 
nately, it is difficult to scale this to optical frequencies, 
where magneto-optic effects are weak. Several related 
systems have subsequently been studied [6UTT]. Of par- 
ticular interest is the system proposed by Hafezi et al. [7] , 
which consists of optical ring resonators arranged in a 
lattice, like a 2D version of the CROW (coupled res- 
onator optical waveguide) [12]. Additional waveguides 
are used to couple the resonators, so that each resonator 
mode couples only (or mainly) to neighboring resonator 
modes with the same direction of rotation. In the tight- 
binding (weak coupling) regime, phase shifts in the inter- 
resonator couplings were shown to play the role of a mag- 
netic vector potential in the tight-binding hopping am- 
plitude, carrying opposite signs for the two directions of 
rotation, which play the role of spin. The couplings may 
be tuned so that the "magnetic field" is uniform, where- 
upon each spin sector exhibits Landau levels in a Hofs- 
tadter butterfly spectrum [13] , with one-way edge modes 
which are topologically protected against spin-conserving 
perturbations. Although the system is reciprocal (time- 
reversal maps the two spin sectors onto each other), and 
therefore cannot be used for optical isolation, Hafezi et 
al. have suggested that its topologically-protected edge 
states may be used as robust optical delay lines [7]. 

In this Letter, we discuss some subtle and interesting 
features of the system of Hafezi et al. which do not appear 
in the tight-binding analysis. In the "zero field" regime, 
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where the synthetic magnetic vector potential induced 
by the inter-resonator couplings vanishes, there exists a 
topological insulator phase at large coupling strengths. 
In this phase, each spin sector exhibits an analog of the 
zero- field quantum Hall effect, including one-way edge 
states, much like the prototypical 2D topological insu- 
lator model of Kane and Mele [14]. Adding an integer 
number of magnetic flux quanta per unit cell maps the 
band structure back onto the zero field case. Our results 
are obtained using the transfer matrix method, which has 
previously been applied to the CROW [I5j [16] , and has a 
much wider domain of validity for such systems than the 
tight-binding method [16 . These calculations are also of 
interest in that they yield topological insulator physics 
without the use of tight-binding, or even of a Hamilto- 
nian. Finally, as a straightfoward extension of the trans- 
fer matrix analysis, we study the band structure of the 
PT-symmetric (parity /time-reversal symmetric [17]) lat- 
tice, which may be realized by adding equal and opposite 
amounts of gain and loss to parity-symmetric regions of 
the coupling waveguides. PT-symmetric photonic lat- 
tices are known to have several unusual properties, such 
as bifurcations between real and complex bands [T8H2T] : 
they have recently been realized experimentally using op- 
tical fiber networks [22] , and the present system could be 
realizable in a similar way. We show that by adding PT- 
symmetric gain and loss to the couplers, we can achieve 
directional amplification, amplifying the one-way state 
on one edge and damping the counter-propagating state 
on the opposite edge. 

We consider the square lattice shown in Fig. [I] At each 
lattice site is a ring resonator, whose modes are coupled 
to the modes on neighboring lattice sites having the same 
"spin" (direction of rotation within the resonator). As 
proposed in Ref. [7 and depicted in Fig. [TJc), this cou- 
pling can be implemented by using a waveguide to join 
the rings. However, for our present purposes it will be 
useful to employ a more abstract representation for the 
couplings. We take the spins to be completely decou- 
pled, and restrict our attention to a single spin sector. 
Let n = (x n ,y n ) denote a lattice site, n + x denote the 
site one unit in the +x direction, n + y the site one unit 
in the +y direction, etc. We specify the coupling be- 
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FIG. 1. Coupled optical resonators in a square lattice, (a) 
Schematic showing the couplings between neighboring loop 
resonators, including the wave amplitudes which enter into 
the coupling relations |l]). (b) Schematic of the lattice over 
several periods, (c) A typical waveguide loop which can give 
rise to the couplings shown in (a). 



tween the resonator at n and the one at n 
complex coupling coefficients r nxi r' r 
ilarly, we denote the coupling between n and n + y by 
another four coupling coefficients r nyi r' ny , t ny , and t' ny . 
These relate the wave amplitudes in the resonator — see 
Fig. [IJa) — according to 
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Here, and in the following, the dummy index \i may take 
the value x or y. The parameter <p is the phase delay 
across each quarter of the resonator (assumed equal). 
The matrices S n/Jj have the form of scattering matrices, 
with {r nfI ,r f nfI } and {t nfl ,t f nfl } playing the role of "reflec- 
tion" and "transmission" coefficients respectively. They 
express the most general possible form of linear spin- 
conserving coupling between the resonators. In princi- 
ple, each of the coupling coefficients {r n/x , r^ M , t n/x , 
can be independently varied by appropriately tuning the 
underlying waveguide coupling parameters and phase de- 
lays [23] . In an experimental system, both (j) and the cou- 
pling coefficients will depend implicitly on the operating 
frequency. However, here we will treat them as indepen- 
dent quantities; when calculating the bandstructure, (j) 
will play the role of frequency (or energy) [IS] . 

Now consider the special case where the "reflection" 
coefficients r n/i and r f nfl vary between different cells by 
some phase factor, according to 
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Here, and A 1 ^ are functions of the lattice coordinates 
n, playing the role of a synthetic magnetic vector poten- 
tial. These gauge relations generalize the criterion used 
in Ref. [7], which involved phase differences in the tight- 
binding hopping amplitudes. 

Suppose the vector potential is associated with a 
uniform rational magnetic flux through each unit cell: 
A n + AV n+x ~ K+ y ~ A l = 2ttP/Q, where P and Q are 
integers. When Q = 1, i.e. for any integer number of flux 
quanta through each unit cell, we can show that the band 
structure is the same as in the "zero field" (A* = A 1 ^ = 0) 
system. In this case, the magnetic unit cell has the same 
periodicity as the lattice, and the system supports solu- 
tions of the form 

a n+ , = e^ +A ^a n , b n+ , = e^+ A ^b n . (4) 

Inserting (|3| into ([T]) yields the quadratic equation [23] 

e -4# _ Be -2%4> _ c = 0, (5) 
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As we shall see, for unitary couplings this gives four 
real bands in the periodic space <p G [— 7r,7r]: two in 
[— 7r/2,7r/2] obtained by directly solving (|5|, and the 
other two by adding ±7r. Since the vector potential does 
not appear in the dispersion relation, the band structure 
for any integer number of flux quanta per unit cell is the 
same as for the zero-field system. This property relies 
crucially on the fact that in ^ there is no phase vari- 
ation in the {t nfI ,t f nfI } coefficients. For the waveguide 
coupler shown in Fig. [ljc), this condition is satisfied if 
the sum of the phase delays on the two arms of the cou- 
pler is kept constant [23] . 

The behavior of the zero- field (or integer flux) system is 
not as simple as one might initially suppose. The system 
possesses an inherent degree of non-reciprocity in each 
spin sector, arising from the coupling equations 0; as it 
turns out, this gives rise to a topological insulator state 
(or, within each spin sector, a zero-field quantum Hall ef- 
fect [3]). The rest of this Letter will focus on this case. As 
for non-integer fluxes (Q ^ 1), the present approach gives 
essentially the same physics as the tight-binding analysis 
of Ref. [7J: using the Landau gauge A^ = (P/Q)y n and 
A v n = 0, we can define a Q x 1 magnetic unit cell in which 
the amplitudes obey 
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and similarly for b. This gives rise to AQ bands, analo- 
gous to Landau levels in the two-dimensional magnetized 
electron gas. The band structure may be obtained by a 
multi-cell calculation, and a Hofstadter butterfly spec- 
trum is obtained by continuously varying the magnetic 
flux [71CE3]. 
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If the couplings conserve energy, the coupling matrices 
are unitary: = S^ 1 . In this situation, we expect the 
band structure 4>(K X , K y ) to be real (for K x and K y real). 
This is easily proven by using the parameterization 
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where M G [0, 7r/2] and X/i? W £ [0? 27r]. For the zero- 
field system, Eq. (pi then simplifies to 
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For real 

|y| <sin(^ + ^) < 1, 
and hence Eq. ([9| yields real bands of the form 
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The above calculation also yields the phase diagram. 
Band-crossing points occur if, and only if, the inequality 
(10) can saturate, i.e. 6 X + y = 7r/2, or equivalently 
rx\ + \ r y\ 2 = 1- This defines a boundary between two 
distinct insulator phases, one of which turns out to be 
topologically non-trivial. In order to demonstrate this, 
we specialize to cp^ — X\i — 0, £^ = 7r/2, and 6 X = 6 y = 0, 
so that 5^ is symmetric (£ M = t' = icos#), and 
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The projected band diagram for a semi-infinite "lattice 
strip" (which has a finite width in the y direction and is 
infinite in the x direction) can be calculated in a similar 
way [23 . Fig. [2] shows the results for several values of 
0. For < 7r/4, the band gap contains edge states with 
both group velocities on each edge, corresponding to a 
trivial insulator. For 6 > 7r/4, the system is a topological 
insulator: within the given spin sector (clockwise), there 
is one single positive group velocity mode confined to 
the upper edge, and one negative group velocity mode 
confined to the lower edge (Fig. |3|. 

We note that these results would not show up in the 
tight-binding analysis of Ref. [7 , where tight-binding 
amplitudes were assigned to the resonator modes. The 
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FIG. 2. Projected band diagram for the semi-infinite res- 
onator lattice, with 20 lattice sites in the y direction. The 
coupling coefficients are given by Eq. (J7| with ip^ — x^ — 0, 
£n = tt/2, and X — y — 0. The synthetic gauge field defined 
m Eq. Q is zero. The bands are degenerate when O — n/A 
(upper right-hand figure). For > 7r/4, each band gap is 
populated by a pair of one-way edge states. In the = OAtv 
plot, the points labeled A and B, located at <j) — 7r/4, indicate 
the edge states whose intensity profiles are shown in Fig. [3] 




FIG. 3. Variation of edge state intensity with y lattice coordi- 
nate, demonstrating single-edge confinement and exponential 
decay into the bulk. The two edge states, labeled A (filled 
circles) and B (open circles), have equal cj) — tt/4, and occur 
at K x — =p 1.587 respectively, with opposite group velocities. 
The lattice parameters are those used in Fig. [2] with = OAtv. 
We are considering the clockwise spin sector, as depicted in 
Fig. [I] The intensities are defined as the mean of the squared 
wave amplitudes on the four arms of the resonator at each 



lattice site, i.e. 



+ K\ 2 + \c n \ 2 + \d n \ 2 )/4. 



topological insulator phase, 6 X + y > 7r/2, corresponds 
to strong coupling between the resonators, rather than 
weak coupling. Furthermore, the one-way edge states 
arise from the internal structure of the modes at each 
lattice site. Intuitively, each edge state corresponds to a 
mode which undergoes critical coupling from resonator 
to resonator along one edge, which only occurs if the 
propagation direction along that edge coincides with the 
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• Re[jRf*] < 0, lower edge 
— Re[Ka] > 0, upper edge 
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FIG. 4. Amplification and damping of the edge states in- 
duced by the PT symmetric gain/loss parameter T, defined 
in Eq. ( |15[ ). All other coupling parameters are the same as in 
Fig. [2] with = 0.47T. Both edge states acquire an imaginary 
component — iT to the wave vector, resulting in one edge state 
being damped and the other amplified. Inset: the intensity 
profiles of the two edge states at Y = 0.2, showing that they 
remain well-confined. 



chosen direction of rotation within each resonator. 

Finally, we study what happens when the resonator 
lattice is PT-symmetric. PT-symmetric topological in- 
sulators are probably impossible to realize in condensed- 
matter systems, but may be possible in optical resonator 
systems, where gain and loss can be selectively added to 
the resonators and/or coupling waveguides. The transfer 
matrix method is well-suited for studing this: we sim- 
ply replace the unitary coupling matrices with matrices 
obeying the PT-symmetry relation for scattering matri- 
ces 24] [26]: 

VTS fI VT=S~ 1 . (13) 

Here, V is a parity operator and T is a time-reversal 
operator. We choose V = [0,1; 1,0] and T to be the 
complex conjugation operator. For the setup shown in 
Fig. [TJc), setting ^2 = ipl causes to satisfy (13) with 
this representation of V and T. The coupling coefficients 
can in general be parameterized by 
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with \rr'\ + \t\ 2 = 1. Note that this parameterization 
differs slightly from that of Ref. [26] , which considers re- 
ciprocal systems, for which the scattering matrices are 
symmetric; our S^s need not be symmetric. For con- 
venience, however, we set cp = and cp' = tt/2. The 
coefficients then have the form 



e r sin 



1 cost 



0, t' 



l cos 
-r 



suit 



(15) 



The gain/loss parameter T characterizes the strength of 
the unitarity breaking, with T = corresponding to the 
unitary coupling considered previously. 



Fig. [4] shows the effect of PT-symmetric gain/loss on 
the edge states of the semi-infinite strip. Here, we in- 
troduce non-zero T to the couplings in the x direction, 
while keeping the y couplings unitary [27] . The two edge 
states (which are localized to opposite edges and have 
opposite group velocities) acquire the same imaginary 
component to their wave vectors. Thus, one of them 
is damped and the other amplified. This has a simple 
interpretation in terms of the coupling scheme shown 
in Fig. [ljc): if we add gain to (say) the lower arm of 
the coupling waveguide and loss to the upper arm, the 
backward-propagating edge state (on the lower edge) is 
concentrated in the upper coupler arms, and is damped; 
while the forward-propagating edge state (on the upper 
edge) is concentrated in the lower coupler arms, and is 
amplified. Previous authors have studied PT-symmetric 
waveguides, and demonstrated that modes concentrated 
at different transverse positions can be selectively ampli- 
fied and damped [20]. Such waveguides, however, are re- 
ciprocal, so each amplified (damped) mode has a counter- 
propagating partner which is amplified (damped) by an 
equal amount. In the present system, topological insu- 
lator edge states have no counter-propagating partners 
in the same spin sector. Hence, if we consider the entire 
strip (rather than a single edge) as an optical waveguide, 
the introduction of PT-symmetric gain/loss gives us a di- 
rectional amplifier, or optical diode. Furthermore, unlike 
the PT-symmetric diode studied in Ref. [28], this device 
does not require optical nonlinearity. 
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5 



Supplementary material 



Appendix A: Inter-loop Coupling Coefficients 

In Eq. 0, we abstract away the couplings between resonators into a set of "reflection coefficients" and "transmission 
coefficients" , reproduced here for convenience: 
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In this section, we discuss how these coefficients can be related to the parameters of the underlying waveguide-coupling 
mechanism, such as that discussed in Ref. [7 . As shown in Fig. []Jc), the resonator amplitudes a, 6, c, and d (we 
will omit the redundant subscripts for simplicity) are coupled to waveguide amplitudes a, /3, 7, and 5, by coupling 
relations H6l: 
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Without loss of generality, r, r', k, and k' may be taken to be real. If the phase delays along the two arms of the 
waveguide are ipi and ip2 respectively, then 



a = e^ 2 5, 



p = e^ 1 7. 



(A6) 



By combining Eqs. (A4)-(A6), we obtain 
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Comparison with ( Al )-( [A3| ) immediately yields the corresponding values of r, r', t, and t f . 

As indicated in Ref. [7 , the synthetic magnetic vector potential can be implemented by altering the phase delay 
ipi on different lattice points, wh ile k eeping ipi + ^ 2 fixed. This leaves the t and t' coefficients, which are given by the 
off-diagonal matrix elements in (A7), independent of ^1, in agreement with Eq. (3). The relevant component of the 
resulting vector potential is A = 



Appendix B: Gauge Transformations 



In this section, we provide additional details a bou t the gauge structure of the resonator lattice. The coupling 
relations between resonator amplitudes, given by (A1)-(A3), can be simplified by eliminating the c and d variables, 
which yields 
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When the coupling coefficients obey the gauge relations (3), the system of equations (Bl) simplifies to 
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We look for solutions to these two equations, using the Bloch ansatz Q. When writing down the amplitudes for 
negative displacements and/or composed displacements using this ansatz, some care is needed to ensure that the A's 
are evaluated at lattice positions consistent with Q. For example, 
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Plugging these relations into ( B2 ) yields 
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The resulting dispersion relation (0 versus K x and K y ) is independent of and A^. 
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Appendix C: Projected band diagram 

The calculation of the projected band diagram for the resonator lattice is much like the calculation of the CROW's 
dispersion relation by Yariv et al. [16] . Suppose we have a semi-infinite strip, whose width is N lattice sites in the y 
direction and infinite in the x direction, with the coupling matrices independent of the x coordinate. We seek solutions 
of the form 

[a, 6, c, d) n+x = e lKx [a, 6, c, d) n . 

Let us consider a sing le column of the lattice, and replace the site index n with j G {1,2, •• • , N} denoting the y 
coordinate. From (Al), 
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From (A2), 
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for j = 1, • • • , N — 1. At the edges of the lattice, we have the boundary conditions (refer to Fig. 1 of the paper): 
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Combining (C1)-(C3) yields the eigenvalue equation 
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where Ma and Mb are the following 2iV x 2N matrices: 
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To generate the projected band diagram, we sweep through th e de sired range (usually [— 7r/2,7r/2] due to the tt- 
periodicity of the band structure), solving the eigenvalue problem (C4) at each <j). The eigenvalues which do not have 
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unit modulus, which correspond to evanescent modes, are discarded; the arguments of the rest are the value(s) of K c 
for that value of 0. 
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